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ABSTRACT 

We present an X-ray and radio study of the wind nebula surrounding the central pulsar 
> ■ PSR J1811 — 1925 in the supernova remnant Gil. 2— 0.3. Using high resolution data obtained 

with the Chandra X-ray observatory and with the VLA radio telescope we show the X-ray and 
radio emission is asymmetric around the pulsar, despite the latter's central position in the very 
circular shell. The new X-ray data allows us to separate the synchrotron emission of the pulsar 
wind nebula from the surrounding thermal emission and that from the pulsar itself. Based 
on X-ray data from two epochs, we observe temporal variation of the location of X-ray hot 
spots near the pulsar, indicating relativistic motion. We compare thermal emission observed 
within the shell, which may be associated with the forward shock of the pulsar wind nebula, to 
thermal emission from a nearby portion of the remnant shell, the temperature of which implies 
an expansion velocity consistent with the identification of the remnant with the historical event 
of 386 A.D. The measured X-ray and radio spectral indices of the nebula synchrotron emission 
are found to be consistent with a single synchrotron cooling break. The magnetic field implied 
by the break frequency is anomalously large given the apparent size and age of the nebula if a 
spherical morphology is assumed, but is consistent with a bipolar morphology. 

Subject headings: pulsars: general — pulsars: individual (AX J1811.5— 1926, — supernovae: individual 
(Gil. 2-0.3) — X-rays: general 
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1. Introduction 

The overall structure and evolution of the X- 
ray emission produced by a spherically symmet- 
ric relativistic particle wind coming from a young 
pulsar in a supernova remnant has been modelled 
extensively (Reynolds & Chevalier 1984; Kennel 
& Coroniti 1984; van dcr Swaluw et al. 2001; 
Blondin, Chevalier, & Frierson 2001). In addition, 
some progress has been made on the fine structure 
and time dependence of the sites of shock acceler- 
ation in the inner regions (Gallant & Arons 1994). 
The observational tests of these models have been 
dominated by discussions of the Crab nebula. The 
advantages of studying the Crab, with its high lu- 
minosity, known age, and well observed pulsar, are 
many. However, the Crab is less than an ideal 
test case for understanding pulsar wind nebulae 
(PWNe) in supernova remnants (SNRs) as a class 
due to its anomalously high luminosity and lack 
of an observable SNR shell. There are a number 
of observed radio composite SNRs which consist 
of a shell and a spectrally distinct inner nebula, 
presumably PWNe. Only in a few cases, however, 
has the central pulsar been detected (eg. W44, see 
Kaspi & Hclfand 2002 for a review, also Gotthelf 
2001), leaving the spin-down energy driving the 
emission to be inferred from the PWN itself. 

Gil. 2— 0.3 is a remarkably spherical young 
SNR plausibly associated with the historical event 
of 386 A.D. recorded by Chinese astrologers (Clark 
& Stephenson 1977). A 65 ms pulsar with spin- 
down energy E = 6.4 x 10 36 ergss _1 was dis- 
covered in X-rays by ASCA (Torii et al. 1997, 
1999). The characteristic spin-down age is much 
greater than the apparent age of the SNR and 
that inferred from its remarkably central position 
within the remnant given any reasonable trans- 
verse velocity (Kaspi et al. 2001, hereafter Paper 
I). The simplest explanation for this apparent age 
discrepancy is that the initial spin period of the 
pulsar is very near the current spin period, which 
in turn suggests that E has remained nearly con- 
stant since the supernova explosion. The distance 
to Gil. 2— 0.3 has been reasonably well determined 
to be d ~ 5 kpc from HI measurements (Becker, 
Markert, & Donahue 1985; Green et al. 1988). 
We therefore have a system where the assump- 
tions of spherical symmetry, at least in regards to 
the shell around the pulsar, and constant known 



energy input are observationally supported, and 
reasonable estimates of the physical size of various 
components can be made. 

The bright spherical radio shell, reminiscent 
of the remnants of Tycho's and Kepler's super- 
novae (Downes 1984), and the lack of an obvious 
hard central component in the ROSAT PSPC im- 
age, initially led some authors to conclude that 
Gil. 2— 0.3 was the remnant of a Type la super- 
nova (Reynolds et al. 1994). The first hints of a 
central plerionic component came from high fre- 
quency single dish radio measurements which sug- 
gested a central flattening of the radio spectrum 
(Morsi & Reich 1987). ASCA observations first 
clearly demonstrated that there was a central, 
non-thermal X-ray source (Vasisht et al. 1996). 
This component was resolved by Chandra, the 
image of which has been presented in Paper I. 
Further single dish radio observations have clearly 
demonstrated the existence of a flat central radio 
component (Kothes & Reich 2001). Archival VLA 
observations separate the central plerion from the 
shell and measure the spectral index to be a r ~ 
0.25 (defined as 5* = S r (v/v r )~ ar , where S r is 
the energy flux density at frequency v r ), a fairly 
typical value (Tarn, Roberts, & Kaspi 2002, here- 
after Paper II). In this paper we combine the ra- 
dio data with the Chandra data to make spectral 
measurements of the pulsar itself, the hard pul- 
sar wind nebula, and soft thermal emission pos- 
sibly related to the PWN which we compare to 
a portion of the shell. We also present evidence 
for apparent temporal variation in the positions of 
emission features in the central portion of the X- 
ray PWN, and make a detailed comparison of the 
X-ray morphology with a new 3.5 cm VLA image. 
We then compare these observations to models of 
PWN evolution and explore the possibility that 
within the SNR shell a wind termination shock 
produces the variable hard emission, and the for- 
ward shock of the PWN expanding into the stellar 
ejecta produces the interior thermal structures. 

2. Observations and Analysis 

2.1. X-Ray and Radio Observations 

NASA's Chandra X-ray Observatory observed 
Gil. 2— 0.3 at two epochs, the first (Sequence 
Number 50076) on 2000 August 6, and the sec- 
ond (Sequence Number 50077) on 2000 October 
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15, for 20 ks and 15 ks respectively. The remnant 
was positioned on the back-illuminated CCD chip 
S3 of the ACIS instrument in standard exposure 
mode with a time resolution of 3.2 s, too coarse 
to detect pulsations from the pulsar. The data 
were reduced and analysed with CIAO 2.2.1 using 
Caldb 2.12, XSPEC V.ll.1.0 and the MIRIAD 
software packages. Details of the image analysis 
and a three color X-ray image of the SNR can be 
found in Paper I. Here, we focus on the emission 
interior to the shell. 

Analysis of archival 20 cm and 6 cm VLA data 
presented in Paper II revealed the extent of the 
radio PWN by showing that most of the central 
enhancement has a significantly flatter spectrum 
than the surrounding shell. In order to further en- 
hance the relative brightness of the PWN over the 
shell, we have obtained 3.5 cm data with the VLA 
in the DnC, C, and CnB configurations as part 
of a campaign to better constrain the broad-band 
radio spectrum and measure the shell expansion. 
Details of the radio observations will be presented 
elsewhere (Tarn et al. in preparation). 

2.2. Spatial Analysis 

The multi- wavelength image of Gil. 2— 0.3 is 
presented in Figure 1. The 3.5 cm radio image is 
shown in green, the soft (0.6-1.65 keV) X-ray im- 
age, dominated by thermal emission, is in red, and 
the hard (4-9 keV) X-ray image is in blue. This 
image shows the relationship between the interior 
thermal X-ray emission, the radio PWN and the 
bright, hard X-ray emission of the PWN. It also 
shows how the relative brightness of the soft X-ray 
emission and the radio emission changes around 
the shell. The point source at the center is the 
pulsar, which is only seen in the X-ray images. Pa- 
per I noted that the hard X-ray PWN consists of 
a bright region southwest of the pulsar, a narrow, 
almost linear component to the northeast, and a 
larger, faint, diffuse component with reasonable 
bi-lateral symmetry around the pulsar (see Fig- 
ure 2 in Paper I). 

In Figure 2 we show the exposure corrected, 4-9 
keV summed image of the two observing epochs, 
smoothed with a 5" Gaussian, with 6 cm radio 
(which better shows a combination of flatter and 
steeper features) contours overlain. The SNR shell 
predominantly emits soft, thermal X-ray emission 
which is mostly invisible above 4 keV, and so the 



features in this image show regions of non-thermal 
emission. The central PWN is quite hard, and 
the greyscale levels in this image have been chosen 
to highlight the difference between the bright and 
faint parts of the PWN. The extent and shape of 
the radio PWN correlate well with the faint hard 
PWN X-ray component. To the northeast of the 
pulsar, the narrow bright hard X-ray feature goes 
through a region of low radio luminosity, and then 
terminates at a bright arc in the radio emission. 
To the southwest, the bright hard X-rays coincide 
with the radio emission, with a radio enhancement 
coincident with the X-ray peak (the 'spot' region 
referred to below). The shell is clearly seen in 
hard X-rays, although they are concentrated to- 
wards the inner edge of the shell, especially in the 
bright southwestern section. These hard X-rays 
are largely interior to the peak emission of the ra- 
dio shell in this region. 

In Figure 3 we show the soft (0.6—1.65 keV) 
X-ray emission with the same radio contours as in 
Figure 2. There is a large, soft X-ray structure 
which we refer to as the soft "PWN", which is 
centered on the pulsar, reaching to the northwest 
and southeast. The ratio of medium to soft X- 
ray flux, suggests this structure does not have any 
large spectral variations (see Paper I). This soft 
X-ray emission seems to outline much of the ra- 
dio PWN and largely coincides with radio features 
that were shown in Paper II to have steeper spec- 
tra consistent with the radio shell emission rather 
than the PWN emission. The soft X-ray and ra- 
dio shells are quite similar, except the radio emis- 
sion is slightly more extended than are the X-rays 
around most of the shell. However, the bright en- 
hancement in the southeastern portion of the shell 
seems to be broader in thermal X-rays relative to 
the radio than elsewhere in the shell. 

To study the details of the PWN, we produced 
separate hard X-ray images for the two observa- 
tion epochs whose coordinates we aligned by us- 
ing the bright emission of the pulsar to correct 
for slight differences in the aspect solution. These 
were then exposure corrected and smoothed with 
a 3" Gaussian (Figure 4). For comparison, we also 
present the same region from the 3.5 cm radio ob- 
servation, which best shows the PWN due to its 
flatter spectra relative to the shell. The brightest 
features in both X-ray images other than the pul- 
sar are two spots to the southwest of the pulsar, 
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which arc ~ 10a and ~ la above the surround- 
ing PWN level. To better see the relative ampli- 
tudes, we also show line intensity profiles of all 
three images of the slice indicated in the images 
(Figure 5). The positions of the two spots were de- 
termined for each epoch both by centroiding and 
by fitting a 2-D Gaussian on top of a constant 
background, the two methods giving consistent re- 
sults. The latter method gives uncertainty esti- 
mates with which we add in quadrature a relative 
aspecting uncertainty of 0.31" using the pulsar as 
our reference (see aspecting calibration memos by 
M. Markevitch 6 ). D uring the first observation, 
the physical distances of the spot centers from the 
pulsar were 0.20 ±0.01 d 5 pc and 0.34 ±0.01 d 5 pc, 
where d§ — d/(5kpc). By the second observa- 
tion the apparent centers had moved significantly, 
- 3.49" ±0.31" and - 1.90" ±0.44" away from the 
pulsar, implying apparent transverse velocities of 
~ lAd 5 c and ~ 0.8ds c. 

2.3. Spectral Analysis 

2.3.1. Methodology 

Chandra's high resolution imaging capabilities 
allow the measurement of point source spectra to 
be made with remarkably low background levels. 
We initially used a 7 pixel radius circular region 
to extract a pulse height spectrum from the pulsar 
with an annular background region from 7 to 10 
pixels. The spectrum showed an enhancement at 
high energies, a clear signature of pilcup. From 
the observed count rate of ~ 8 x 10 _2 ctss _1 , we 
estimate a pile-up fraction ~ 5%, which is enough 
to significantly affect the measurement of the spec- 
tral index. Around half of the counts are contained 
within one pixel at the center of the PSF, which 
is where most of the pile-up occurs. By excluding 
the central pixel, we can minimize the effects of 
pilcup, but must then worry about the variations 
in the PSF as a function of energy which will af- 
fect the spectrum when a restricted extraction re- 
gion is used. Since the current version of CIAO is 
not able to calculate the effective area when only 
a fractional part of the PSF is enclosed, we ex- 
amined the PSF enclosed energy curves and care- 
fully chose annular extraction regions for the on 
source and background spectra whose difference 

'http://asc.harvard.edu/cal/Hrma/ 
hrma/optaxis /platescale /geom.public. html 



would provide a relatively constant enclosed en- 
ergy fraction over most of the Chandra waveband. 
Choosing an annular source extraction region with 
inner and outer radii of 0.3" — 1.0" , and a back- 
ground annular region of 1.0" — 1.38", we obtained 
an enclosed energy fraction of ~ 0.37 which varies 
by only a few percent in the 1-9 keV range. We fit 
the resulting spectrum to an absorbed power-law 
model, the results of which are listed in Tabic 1. 

For extended sources, both the cosmic X-ray 
and particle background are non-negligible, and 
vary with position on the detector. At the mo- 
ment, there is no consensus on the best way of ap- 
proaching this problem. We have therefore used 
the following procedure, based on one developed 
by Roberts, Romani & Kawai (2001) for analysis 
of ASCA GIS data. We assume that the back- 
ground consists of three components: a particle 
component which is not necessarily spatially cor- 
related with the telescope effective area, an ex- 
tragalactic component which suffers from absorp- 
tion by the Milky Way, and a Galactic compo- 
nent. We further assume that the particle back- 
ground is relatively time independent once the 
data are screened for background flaring events 
and that the two cosmic contaminating compo- 
nents uniformly illuminate the telescope over the 
S3 chip field of view. The Chandra science sup- 
port center has made available blank-sky obser- 
vations of high Galactic latitude fields which we 
assume accurately represent the particle and un- 
absorbed extra-Galactic components of the back- 
ground. 

Using the above assumptions, we account for 
the different background contributions in the fol- 
lowing way. We extract spectra from the region of 
interest and a source-free region of the chip. From 
each of these, we subtract spectra extracted from 
blank-sky data sets from the same chip regions in 
order to remove the particle background contri- 
bution. Since the effective exposure varies as a 
function of chip position, we scale the local back- 
ground spectra by the ratio of average exposures 
of the source and background regions using the ex- 
posure maps which were used in producing the ex- 
posure corrected images of Paper I. Although this 
corrects for the gross differences in effective area 
between the source and background extraction re- 
gions, it should be noted that the variation in effec- 
tive area with position is energy dependent, which 
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is not taken into account with the current proce- 
dure. However, this assumption allows us to ac- 
count for the background in a model independent 
way. In practice, for Gil. 2— 0.3 the spectral re- 
sults are not strongly sensitive to the background 
subtraction method used, as the source to back- 
ground ratio is high. To create effective area and 
response matrix files, we make a binned image of 
the source region photons from which we subtract 
an image created from the blank-sky data set, and 
use this as the input weighting map to the CIAO 
mkwarf and mkrmf tools. The resulting source and 
background spectra and calibration files are used 
in XSPEC to fit models to the spectra. Note that, 
given the assumptions described above, all three 
background components are accounted for by this 
procedure. 

2.3.2. Results of Model Fitting 

In Figure 6, we show the extraction regions 
used for spectral analysis. The greyscale again 
shows the hard X-ray emission, while the contours 
show the soft X-ray emission. The interior of the 
shell has enhanced emission over the general back- 
ground, so for features within the shell, we choose 
a background region interior to the shell. In the 
case of the shell itself (region 1), a region exterior 
to the shell is chosen as background. 

We analyzed the bright region of the shell (re- 
gion 1) near that of a portion of the interior soft 
feature (region 4). In order to compare the two 
regions, we initially used the constant tempera- 
ture, plane-parallel shock plasma model PSHOCK 
in XSPEC for both. However, there is a clear 
hard excess in region 1 seen both in the spectrum 
and the hard X-ray image (Figure 6), while there 
is very little emission above 4 keV in region 4. 
This hard X-ray emission is expected from shock 
acceleration of particles in SNR shells, the pre- 
sumed source of Galactic cosmic rays. To account 
for this, we add the synchrotron roll-off SRCUT 
model (Reynolds & Keohane 1999), which simu- 
lates the roll-off of the synchrotron emission spec- 
trum of an exponentially cut off power-law distri- 
bution of electrons, to the region 1 fit, fixing the 
radio spectral index a — 0.56 as determined in 
Paper II and the 1 GHz flux extrapolated from 
the 1465 MHz image. We then fit for the roll-off 
break frequency which will reproduce the observed 
hard excess. In both regions, the spectral features 



seemed to be shifted by ~ 1.5 channels. Since 
we used a 10 eV/channel binning scheme when 
generating the response matrices, this is similar 
to the +16 eV zero point gain shift reported on 
31 Jan. 2002 to the Chandra users committee 7 . 
We corrected for this by applying a 1.6 channel 
bin shift to the response matrix. A significant im- 
provement to the region 4 fit could be made by 
adding a weak narrow line at 1.18 keV, and there 
was weaker evidence for this feature in the other 
spectra. The strength of this line varied with the 
size of the extraction region, and not the total flux. 
We therefore suspect that this is some background 
component not adequately removed by the proce- 
dure described in the above section. 

For region 4, (and subsequently regions 2 and 
3), we were able to obtain a reasonable fit using the 
simple absorption and shock models, but found 
the much greater signal to noise region 1 spectrum, 
with ~ 25, 000 counts, showed significant residu- 
als around all of the spectral line features. Using 
the VPSHOCK model, which allows the individ- 
ual element abundances to vary, greatly improved 
the fit, but there still seemed to be systematic off- 
sets in the location of the spectral features. Al- 
lowing a rcdshift improved the fit statistically, but 
the implied velocities were too large to be plausi- 
ble. We then tried allowing variable abundances in 
the absorption using the VPHABS model, which 
significantly improved the overall fit resulting in 
a reduced x 2 ~ 1. However, there is significant 
correlation between the individual element abun- 
dances as well as between the abundances in the 
shock model and their corresponding absorption 
abundances. It is also true that calibration un- 
certainties in the response, which are known to 
lead to overestimates of line widths, will have had 
an unknown effect on the abundance values. We 
therefore do not report the individual abundance 
values quantitatively and instead make only gen- 
eral comments. It should be noted that the vari- 
ous model changes only slightly alter the fit tem- 
perature and overall absorption column of the re- 
gion 1 spectrum. The VPSHOCK and VPHABS 
also improved the region 4 spectral fit, but the 
freedom allowed by poorer statistics (from only 
~ 4000 counts in the spectrum) could not effec- 
tively constrain the parameters, and so we only 

7 http://asc. harvard.edu/udocs/ucom.jan02. html 
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report the results of the fits without individually 
varying abundances. 

The resulting temperatures of regions 1 and 4 
were similar, with the shell having abundances 
several times solar. However, these are all sensitive 
to the fit absorption abundances, which tended to 
be high in Si and S, but low in Mg. High resolution 
spectra very sensitive to low energies would be re- 
quired to disentangle the absorption and emission 
features. There were several local minima in the 
multi-parameter \ 2 space near the best fit values, 
with moderately different temperatures (ranging 
from kT = 0.5 — 0.6 keV) and absorptions. The 
ionization timescale upper limit t u parameter of 
the shock model and the SRCUT break frequency 
were somewhat more sensitive, but did not vary by 
more than a factor of 2. We emphasize that the 
90% multi-parameter confidence regions presented 
in Table 1 are statistical around our best fit values 
and may not represent the true uncertainty. How- 
ever, the systematic uncertainties are not likely to 
affect the qualitative interpretation of the results 
presented in the discussion below. 

We also tried fitting the shell portion using 
the VNPSHOCK model which allows for differ- 
ent ion and electron temperatures. There were no 
significant differences between the two tempera- 
tures, indicating that the ions and electrons are 
in rough equilibrium. However, the dependence of 
the model on changes in ion temperature is weak 
(Borkowski, Lyerly, & Reynolds 2001), and mod- 
erate temperature differences could not be ruled 
out. 

In order to determine the synchrotron emission 
properties of the hard PWN (regions 2 and 3), 
the soft thermal emission in the central regions 
near the pulsar needs to be accounted for. We use 
our best fit PSHOCK and PHABS values from the 
interior soft emission (region 4) spectral fits, and 
allow only the normalization and total ur to vary. 
To this we add a simple power-law model to fit the 
hard PWN components. The synchrotron X-rays 
are clearly seen in the data (Figure 7). The results 
of all the fits are listed in Table 1 with 90% multi- 
parameter confidence regions in parentheses. The 
observed fluxes with la errors are listed for the 
energy ranges indicated for the pulsar and regions 
1 and 4. For regions 2 and 3, the unabsorbed fluxes 
of the power-law model component only are given. 



3. Discussion 

Models of PWN evolution in SNRs suggest 
there are two major epochs in the life of the sys- 
tem. In the early evolution, the PWN expands 
into the freely expanding ejecta of the SNR. The 
highly relativistic pulsar wind is terminated at a 
shock near the pulsar, which continuously injects 
electrons and ions having a power-law distribu- 
tion of energies, as well as magnetic field, into 
a surrounding cavity. The particles then gyrate 
in the nebular magnetic field, causing the syn- 
chrotron emission. This cavity itself expands at 
supersonic speeds relative to the surrounding SNR 
ejecta. Therefore, there should be a forward shock 
at the edge of this cavity, similar to the SNR for- 
ward shock moving into the surrounding ISM. 

As the SNR shell expands, sweeping up more 
of the interstellar medium, a reverse shock is 
launched which travels towards the SNR center, 
carrying information about the external medium 
to the inner portions of the remnant. When the 
reverse shock encounters the forward shock of the 
PWN, the PWN is compressed, and reverberates 
several times (van der Swaluw et al. 2001). The 
compressed PWN finally settles into a smaller size 
relative to the SNR shell, and begins to expand 
again, this time sub-sonically. This corresponds 
to the SNR entering the Sedov expansion phase, 
where the dynamics are dominated by the swept- 
up ISM mass. 

Hydrodynamical models of spherically symmet- 
ric PWNc in SNRs suggest the collision between 
the reverse shock and the PWN occurs after ~ 
2000 yr (van der Swaluw et al. 2001; Blondin, 
Chevalier, & Frierson 2001). Although the energy 
input of the pulsar in these models was appropri- 
ate for the high energy output of the Crab and not 
the relatively steady low energy output of PSR 
J1811 — 1925, the time of collision is largely de- 
termined by the propagation of the reverse shock 
which only depends on the dynamics of the su- 
pernova ejecta. We therefore expect the timescale 
in the case of Gil. 2— 0.3 to be roughly similar. 
If we assume the age of the SNR to be 1615 yr, 
corresponding to the 386 A.D. historical associa- 
tion, then we would expect the time of collision 
to be near. If the collision has not yet occurred, 
the PWN should be near its largest extent relative 
to the SNR shell, and its morphology should be 
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largely determined by the geometry of the pulsar 
wind. If the collision has occured or is currently 
occurring, then the PWN should be compressed 
and distorted by the reverse shock. 

The separation of the PWN from the surround- 
ing shell in both the X-ray and radio images sug- 
gests that the reverse shock has not yet reached 
the PWN. The hard X-ray emission concentrated 
at the inner edge of the shell may result from the 
recent passage of the reverse shock, suggesting the 
latter is still near the observed shell. This hypoth- 
esis is somewhat supported by the rather large size 
of the radio PWN (Paper II). In this case, one 
might hope to see evidence of the PWN forward 
shock in the form of thermal X-rays. The ther- 
mal emission interior to the shell, which we refer 
to as the soft "PWN", may be a result of this 
shock. The smaller ionization timescale of this 
feature compared to the nearby shell emission may 
support this interpretation. The similar temper- 
atures suggest the shock speeds are nearly iden- 
tical, which might be unexpected if the radius of 
the soft "PWN" is only 2/3 the radius of the shell, 
as it appears in the image. However, the SNR for- 
ward shock will have been significantly slowed by 
the accumulation of ISM material over the life- 
time of the remnant, and the relative magnitude 
of the velocity jumps at the two shock fronts is 
model dependent (van der Swaluw et al. 2001). It 
is therefore difficult to say what relative shock ve- 
locities, and hence relative temperature, should be 
expected. 

Following the reasoning of Reynolds et al. 
(1994), the Rankinc-Hugoniot shock speed can 
be inferred from the best fit temperature of 
kT ~ 0.58 keV of the shell region to be v s = 
(16kT/3fim p )^ 2 ~ 700 km s" 1 assuming electron- 
ion equilibrium and that the mean mass per par- 
ticle /j,m p is 0.6 times the proton mass. Since the 
ions heat the electrons, the observed electron tem- 
perature is a lower bound on the ion temperature 
and hence the shock velocity. However, the results 
of the VNPSHOCK model suggest that the elec- 
trons and ions are not far from equilibrium. The 
velocity thus derived implies a current expansion 
rate of 9 < 0.03"/<i5yr -1 , giving a free expansion 
age of Tf e — R snr /v s < 5000<i5 yr for the SNR 
radius R 150". If we adopt the Sedov rela- 

tion for a blast wave of radius R$ = 2.5v s t we 
obtain a Sedov age estimate of ts < 2000d 5 yr, 



close to the historical age r = 1615 yr. The ion- 
ization timescale parameter from the VPSHOCK 
model implies a reasonable electron density of 
n es ~ 10cm~ 3 for a 1615 yr old shock, while 
the synchrotron roll-off break frequency implies a 
maximum electron energy of ~ 20 TeV for a typi- 
cal magnetic field strength B ~ 10 [iG (Reynolds 
& Keohane 1999). 

The pulsar has a hard spectrum with no sign 
of a thermal component, although that is not sur- 
prising given the high level of absorption. The 
photon index T ~ 1.3 is similar to what has been 
found recently from other Chandra observations of 
young pulsars (Lu et al. 2002; Pavlov et al. 2001; 
Gaensler et al. 2002). The unabsorbed flux, F x ~ 
4.2 x 10 _12 ergscm~ 2 s _1 in the 1-10 keV energy 
band, is similar to the ~4x 10~ 12 crgscm~ 2 s _1 
quoted by Torii et al. (1997) for the total pulsed 
flux, suggesting that the pulsed fraction is nearly 
100%. 

3.1. Spots as Wisps 

The spatial variability observed in the bright 
PWN spots may be the equivalent of the Crab's 
wisps, which are seen to form at the radius of the 
inner ring (r ~ 0.14 pc) then move outwards at 
a speed v w ~ 0.5c (Mori et al. 2001). In the 
Gallant and Arons (1994) model, the wisp spacing 
roughly corresponds to the post-shock ion Larmor 
radius. In such a scenario, the spots should form 
at the characteristic radii and then propagate out- 
ward and dissipate. The coincident enhanced ra- 
dio emission (Figure 4), with its much longer cool- 
ing timescale, suggests that these spots are per- 
sistently found in this region. We also note that 
there appears to be dynamical activity within the 
narrow feature on other side of the pulsar, but the 
variations in the flux levels are too near the noise 
level to unambiguously identify individual knots 
at both epochs. The bright features may also rep- 
resent a collimated, jet-like outflow. In this case, 
we might expect new knots to be ejected near the 
pulsar and move outwards and diffuse. The appar- 
ent velocity of c or even greater and the brightness 
asymmetry around the pulsar suggest this may be 
the case. 

The asymmetry of the radio emission compared 
with the X-ray emission is puzzling. On one side 
of the pulsar, the emission is coincident, while 
on the other, the radio and X-ray emission arc 
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anti-coincident. Such correlations have been ob- 
served in other PWN, most notably around PSR 
B1509— 58 (Gaensler et al. 2002) whose character- 
istic age is similar to the age of Gil. 2— 0.3. The 
central location and corresponding upper limits on 
the velocity of the pulsar (see Paper I) make inter- 
pretations of such features as a trail marking the 
previous passage of the pulsar through the region 
untenable, and argue against a large asymmetry 
in the birth explosion as a causal factor. 

3.2. Magnetic Field 

In this section, we present a very simplistic 
modification to standard PWN magnetic field esti- 
mates to account for non-sphericity. We can iden- 
tify the location of spot 1, r s \, as being down- 
stream from the pulsar wind termination shock 
r si ft, and use this to estimate the downstream 
magnetic field B t of the PWN. If we assume a ra- 
dial relativistic flow and a postshock velocity of 
c/3, then (eg. Chevalier 2000): 



Bi 



1/2 



( e 1/2 \ 
6 X 10~ 5 rv- G, 



ft(n/r s i)d 5 



(1) 

where e t < 1 is the ratio of magnetic energy to 
total energy in the postshock flow. To allow for a 
non-spherical outflow, we include a beaming frac- 
tion at the termination shock ft < 1. 

The current best estimate of the radio energy 
spectral index is a r ~ 0.25 (defined as S v oc v~ ar ) 
with an acceptable range of a r — 0.15-0.3. The 
corresponding range of the X-ray energy index is 
a x = T - 1 = (0.54-0.90) with a best fit value 
ax = 0.73. The difference between the X-ray 
and radio spectral indices Sa r x ~ 0.5 is consis- 
tent with what is expected for a single synchrotron 
cooling break between the radio and X-ray bands, 
although the uncertainty may be as large as 0.25. 
By assuming Sa rX — 0.5, we can find the break 
frequency from the measured radio index and ra- 
dio and X-ray flux densities: 



S 



xv x 



+0.5 



(2) 



We estimate the flux density of the radio nebula 
at 1.5 GHz to be S r (PWN) ~ 0.4 Jy and of the 



spot region to be S r (spot) ~ 0.04 Jy, although 
these could be off by a factor of two (see Pa- 
per II for a discussion of flux uncertainties). For 
the observed range of values for a ri this results 
in estimated upper limits on the break frequen- 
cies of v b {PWN) < 50 GHz and v b (spoi) < 200 
GHz, with values for the nominal spectral index 
a r = 0.25 of ~ 8 GHz and ~ 25 GHz, respectively. 

If we assume constant energy injection of par- 
ticles over the lifetime of the PWN, as indicated 
by the inferred initial spin being nearly the same 
as the current spin period (Paper I), we can esti- 
mate the average magnetic field from the apparent 
break frequency (eg. Chevalier 2000) from: 



D 



4 x 10" 3 



N 



(^/10GHz) 1 /3(i/l614yr) 2 / 3 



G. (3) 



This is very large compared to the post-shock esti- 
mate above (Eq. 1). Even with the most extreme 
values of flux and spectral index, because of the 
weak dependence on the magnetic field in the 
nebula and in the spot region Bx > 10~ 3 G. If the 
magnetic field is constant throughout the emitting 
volume of the PWN, then the maximum magnetic 
field can be inferred from the energy injected over 
the lifetime of the pulsar assuming negligible ra- 
diative losses (eg. Pacini and Salvati 1973): 



B N < 



6e N Et 



1/2 



■ { ) G 



1/2 



(4) 

where ex is the fraction of the spin-down energy 
which ultimately goes into magnetic field (0.5 for 
equipartition) and fx is the ratio of the emitting 
volume to a sphere of radius rjv ~ 40", the dis- 
tance from the pulsar to the radio arc. In order 
to reconcile Equations 3 and 4, we need a nebu- 
lar filling fraction fx ~ 0.01. This is feasible if 
the emission region is bipolar or toroidal rather 
than spherical, and if the region of high magnetic 
field, where most of the emission takes place, is 
somehow confined to the narrow region of high X- 
ray emission and the surface of the wind bubble. 
The observed PWN structure in the 3.5 cm im- 
age may support this picture (see Figure 4). In 
general, ft < fx, since we expect that the out- 
flow would not become much further collimated 
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after the termination shock, but rather the solid 
angle of the flow would tend to increase with ra- 
dius. This suggests that the pulsar wind is initially 
confined to narrow, polar outflows with opening 
angles 6 < 10° in the case of a bipolar nebula or 
a thin, equatorial sheet in the case of a toroidal 
nebula. 

We introduce the parameters (3 = (Bn / 'B t ) 2 (e t / 'ejv) 
and / = /t//jv < 1- It can be shown that the ar- 
guments of Rees & Gunn (1974) (see also Kennel 
& Coroniti 1984) for the build up of the magnetic 
field to approximate equipartition strength may 
be applied to both a conical and nozzle type flow, 
in which case we expect (3 ~ 1. We can then 
combine Equations 1 and 4 to derive the following 
estimate of the termination shock radius: 



n ~ \0r N ) 



1/2 



20, 



(5) 



suggesting r t should be a few arcseconds from the 
pulsar. In that case, monitoring of the emission 
by Chandra should resolve the region of spot for- 
mation. 

Bandiera, Pacini, and Salvati (1996) noted that 
the extrapolated ASCA PWN spectra combined 
with previous estimates of the radio PWN flux 
led to an anomolously high magnetic field, and 
suggested a period of rapid spin-down early in 
the pulsar's history where the magnetic field was 
dumped into the surrounding region, and confined 
by the external SNR shell. This was previous 
to the discovery of pulsations, and their models 
greatly overestimate the pulsar spin period. It is 
interesting to note, however, that an initial spin 
period Pq ~ 20 ms, similar to that of the Crab 
pulsar's, would provide the necessary energy for 
the magnetic field. Nevertheless, the arguments 
made above for non-sphericity in the nebula obvi- 
ate the need for an ad-hoc invocation of an early, 
short-term, very rapid spin-down era. 

Inferring the magnetic field by extrapolating 
the total PWN radio and X-ray spectra assuming 
a single break may be questionable. There may be 
a change in spectral index over the PWN to which 
we are not sensitive. There is also a possibility 
that the reverse shock has begun to encounter the 
PWN, temporarily enhancing the magnetic field 
at the edges. However, the spot region shows ev- 
idence of current particle injection and is likely 



not yet affected by the reverse shock. The break 
frequency of this region is ~ 5 times higher here 
than for the PWN as a whole, which only lowers 
the magnetic field estimate by a factor of 2, assum- 
ing this feature has existed for the entire lifetime 
of the PWN. 

4. Conclusions 

We have measured the X-ray spectrum of vari- 
ous components of the shell and pulsar wind neb- 
ula of Gil. 2— 0.3 and compared them with the 
radio spectrum and morphology. The shell tem- 
perature and hard X-ray structure are consistent 
with a ~ 2000 yr old remnant nearing the onset 
of the Sedov phase. The pulsar wind nebula itself 
shows evidence of relativistic dynamic evolution of 
bright X-ray enhancements near the pulsar. We 
argue that the reverse shock has not yet reached 
the expanding pulsar wind bubble, and therefore 
the observed structure is largely determined by the 
structure of the wind itself. The morphology and 
magnetic field suggested by the low spectral break 
energy suggests the outflowing wind is highly non- 
spherical. The wind may be largely constrained 
to a narrow, bipolar outflow, or possibly a thin, 
equatorial sheet. Further high-resolution obser- 
vations at X-ray, radio and infrared wavelengths 
could identify the region of spot formation, help 
determine the true morphology, and further con- 
strain the break energy of the spectrum. 
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Fig. 1. — Comparison of X-ray PWN emission and 
radio emission. Red represents photons of energies 
0.6-1.65 keV, green is 3.5 cm radio, and blue is 4-9 
keV photons. 
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Fig. 2. — Comparison of hard X-ray emission and 
radio emission. Greyscale is 4-9 keV X-rays, con- 
tours are 6 cm radio emission. 
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Fig. 3. — Comparison of soft X-ray emission and 
radio emission. Greyscale is 0.6-1.65 keV X-rays, 
contours are 6 cm radio emission. 
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Fig. 4. — ACIS observing epochs 1 (top) and 
2 (bottom), smoothed with a 3" Gaussian, and 
3.5 cm (bottom) VLA X-band image with a 3.1" x 
2.8" restoring beam. The line in the images shows 
the 1-D slice for the profiles shown in Figure 5. 
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Fig. 5. — Intensity profile of a line through the 
spots for the two observing epochs and the radio 
image, showing the displacement between the two 
epochs. The y-axis is in arbitrary units of flux 
density. 
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Fig. 6. — Hard (4-9 keV) X-ray image with soft 
(0.6-1.65 keV) X-ray photon contours, showing 
spectral extraction regions. 
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Fig. 7.— ACIS S3 spectrum of PWN (region 2), 
showing residuals to the absorbed PSHOCK + 
power-law model. The Mg, SI, and S line emis- 
sion from the contaminating thermal emission are 
clearly evident, as is the high energy power-law 
tail from the synchrotron nebula. 
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Table 1 
Spectral fits 



Model Parameter Pulsar Shell (1) PWN (2) Spots (3) Soft "PWN" (4) 



(V)PHABS a 


n H (10 22 cm- 2 ) 


2.36(1.75-3. 19) b 


1.71(1.70-2.36) 


2.14(2.04-2.25) 


2.28(1.86-2.77) 


2.15(2.03-2.27) 


POWER" 


r 


1.11(1.00-1.48) 




1.73(1.54-1.90) 


1.73(1.55-1.95) 




(V)PSHOCK d 


kT (keV) 




0.581(0.553-0.595) 


0.583(frozcn) 


0.583(frozcn) 


0.583(0.517-0.691) 




T U (10 11 scm~ 3 ) 




6.7(5.7-8.1) 


4.2(frozen) 


4.2(frozen) 


4.2(2.3-14.9) 


SRCUT° 


a 

v b (10 7 GHz) 




0.56 
1.80(1.68-1.90) 










Sigh, (Jy) 
ergs cm - 2 s _1 ) 




2.0 








Fx (10- 12 


3.29 ± 0.14 (1-10 keV) 


3.70 ± 0.03 (0.7-8 keV) 


4.44 ± 0.12 f (l-10 keV) 


0.87 ± 0.04 f (l-10 keV) 


0.37 ± 0.01 (0.7-8 keV) 



a Photoclcctric absorption model, allowing the elemental abundances to vary for the shell (1) region only. 
b Values in parentheses are 90% multi-parameter confidence regions. 
c simplc photon power law 

d Plane Parallel shock model of Borkowski, Lyerly, & Reynolds (2001), allowing the elemental abundances to vary for the shell (1) region only. 
G Synchrotron roll-off model of Reynolds & Keohanc (1999). 
f Unabsorbed flux of power law component only. 



